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The nonperturbative approach to soft high-energy hadron-hadron scattering, based on the ana- 
lytic continuation of Euclidean Wilson-loop correlation functions, makes possible the investiga- 
tion of the problem of the asymptotic energy dependence of hadron-hadron total cross sections 
by means of lattice calculations. In this contribution we compare the lattice numerical results to 
analytic results obtained with various nonperturbative techniques. We also discuss the possibility 
to obtain indications of the rise of hadron-hadron total cross sections with energy directly from 
the lattice data. 



The XXIX International Symposium on Lattice Field Theory - Lattice 2011 
July 10-16, 2011 

Squaw Valley, Lake Tahoe, California 



* Speaker. 

"'"Supported by MICINN under the CPAN project CSD2007-00042 from the Consolider-Ingenio2010 program, and 
under the grant FPA2009-09638. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it / 



Hadron-hadron total cross sections and soft high-energy scattering on the lattice 



Matteo Giordano 



1. Introduction 

The problem of predicting total cross sections at high energy from first principles is one of 
the oldest open problems of hadronic physics (see, e.g., [jl|] and references therein), not yet satis- 
factorily solved in QCD. This problem is part of the more general problem of high-energy elastic 
scattering at low transferred momentum, the so-called soft high-energy scattering. As soft high- 
energy processes possess two different energy scales, the total center-of-mass energy squared s 
and the transferred momentum squared t, smaller than the typical energy scale of strong interac- 
tions (|? | < 1 GeV 2 <C s), we cannot fully rely on perturbation theory. A genuine nonperturbative 
approach in the framework of QCD has been proposed in [Q] and further developed in a number 
of papers (see, e.g., for a list of references): using a functional integral approach, high-energy 
hadron-hadron elastic scattering amplitudes are shown to be governed by the correlation function 
(CF) of certain Wilson loops defined in Minkowski space. Moreover, as it has been shown in [^]- 
U, such a CF can be reconstructed by analytic continuation from the CF of two Euclidean Wilson 
loops, that can be calculated using the nonperturbative methods of Euclidean Field Theory. 

In [[7], |8]] we have investigated this problem by means of numerical simulations in Lattice 
Gauge Theory (LGT). Although we cannot obtain an analytic expression in this way, nevertheless 
this is a first-principle approach that provides (within the errors) the true QCD expectation for 
the relevant CF. In this contribution, after a survey of the nonperturbative approach to soft high- 
energy scattering in the case of meson-meson elastic scattering, we will present our numerical 
approach based on LGT, compare the numerical results to the existing analytic models, and discuss 
the possibility to obtain indications of the rise of total cross sections directly from the lattice data. 



2. High-energy meson-meson scattering and Wilson-loop correlation functions 

We sketch here the nonperturbative approach to soft high-energy scattering (see [f7|] for a more 
detailed presentation). The elastic scattering amplitudes of two mesons (taken for simplicity with 
the same mass m) in the soft high-energy regime can be reconstructed, after folding with the appro- 
priate wave functions, from the scattering amplitude ^ud) °f two dipoles of fixed transverse sizes 
^l,2J_> an d fixed longitudinal-momentum fractions f\^ of the two quarks in the two dipoles [0]. 
In turn, the dipole-dipole (dd) scattering amplitude is obtained from the (properly normalised) CF 
of two Wilson loops in the fundamental representation, defined in Minkowski spacetime, running 
along the paths made up of the quark and antiquark classical straight-line trajectories, and thus 
forming a hyperbolic angle % (~ log(.s/ra ) at high energy). The paths are cut at proper times 
±T as an IR regularisation, and closed by straight-line "links" in the transverse plane, in order to 
ensure gauge invariance. Eventually, the limit T — > °° has to be taken. 

It has been shown in [|3|-^] that the relevant Wilson-loop CF can be reconstructed, by means 
of analytic continuation, from the Euclidean CF of two Euclidean Wilson loops, 

M9;T;zx; 1,2) = -4=J* ^f- - 1 , < 2 T) = — tr <^ Texp -ig i A^dx^ \ , (2.1) 

where (. . .)e is the average in the sense of the Euclidean QCD functional integral, and "1 [2]" stand 
for "^i[2]±)/i[2l"- The Euclidean Wilson loops W l 2 are calculated on the following straight-line 
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paths, 

^ : X 1 "® (t) = z + — t + fi®R x , V 2 : X 2 ^ (t) = — t + % 2 , (2.2) 
m m 

with T € [— T,T], and closed by straight-line paths in the transverse plane at T = ±7\ The four- 
vectors pi and p2 are chosen to be p^ 2 ] = fra^dqsin^Oj^cos ^) (taking X4 to be the "Euclidean 
time"), 6 being the angle formed by the two trajectories, i.e., p\ - p% = ra 2 cos 6. Moreover, R[ = 
(0,R i± ,0) and z = (0,z±,0). We define also the CF with the IR cutoff removed as ^ E = lim & E . 
The dd scattering amplitude is then obtained from ^ by means of analytic continuation as 

JZ {dd) {s,t;\,2) = -i2s J d 2 z±e^tf E (e -> -i% ^-ilog(s/m 2 );zr, 1,2) , (2.3) 

where 5 = (p\ + P2) 2 and f = — \q\_\ 2 (q±_ being the transferred momentum) are the usual Man- 
delstam variables (for a detailed discussion on the analytic continuation see where we have 
shown, on nonperturbative grounds, that the required analyticity hypotheses are indeed satisfied). 
In the following, without loss of generality pp, we will take the longitudinal-momentum fractions 
/i,2 = 2> an d suppress the dependence on f\ :2 in and 

3. Wilson-loop correlation functions on the lattice 

The gauge-invariant Wilson-loop CF <S E is a natural candidate for a lattice computation, but 
the explicit breaking of 0(4) invariance on a lattice requires special care. As straight lines on a 
lattice can be either parallel or orthogonal, we are forced to use off-axis Wilson loops to cover a 
significantly large set of angles [|7|]. To stay as close as possible to the continuum case, the loop 
sides are evaluated on the lattice paths that minimise the distance from the continuum paths: this 
can be easily accomplished by means of the well-known Bresenham algorithm JToj] . The relevant 
Wilson loops ;r±;n) are then characterised by the position n of their center and by two 2D 
vectors lu and r±, corresponding respectively to the longitudinal and transverse sides of the loop. 
Setting Wn = W L (fifrij_;d), d = (0,d x ,0), and W L2 = W L (l 2 fr 2J _;0), we define on the lattice 



f L (l lh l 2f X ± )= y* J/ -1, V L (l lh l 2f X x )= Km & L (l n ,l 2 fX ± ). (3.1) 



Here X^ denotes collectively the relevant transverse variables, Xj_ = (d±;ri±,r 2 ±). Also, h\ = 
li\\/Li, where L, = | are defined to be the lengths of the longitudinal sides of the loops in lattice 
units. In the continuum limit, where 0(4) invariance is restored, we expect 

MhwM'M ~ W\T\ = ^,T 2 = %aX x ), V L (h\\,i2\\;X ± ) ^ tf E (6;aX ± ), (3.2) 

where tu\ ■ l 2 \\ = cos 6 defines the relative angle 6 and a is the lattice spacing. 

In [f7[ [Sj] we have performed a Monte Carlo calculation of ^ for several values of the rel- 
ative angle, various lengths and different configurations in the transverse plane. We used 30000 
quenched configurations generated with the SU (3) Wilson action at /3 = 6.0, corresponding to 
a ~ 0.1 fm, on a 16 4 hypercubic lattice with periodic boundary conditions. This choice is made in 
order to stay within the so-called "scaling window": in this sense we are relying in an indirect way 



on the validity of the relation (3.2) between Wilson-loop CFs on the lattice and in the continuum. 
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Figure 1: (Left) The relevant Wilson-loop configuration. Using the 0(4) invariance of the Euclidean theory 
we have put p\ parallel to the x\ axis. (Right) Loop configuration in the transverse plane. 

To keep the corrections due to 0(4) invariance breaking as small as possible, we have kept 
one of the two loops on-axis and we have only tilted the other one as shown in Fig. [I] (left). The 
on-axis loop Wl\ is taken to be parallel to the Xe\ axis, = (Li,0), and of length L\ = 6,8, 
and we have used two sets of off-axis loops Wu_ tilted at cot 6 = 0, ±1, ±2. We have used loops 
with transverse size |ru_| = |?2±| = 1 in lattice units; the loop configurations in the transverse 
plane are those illustrated in Fig. [I] (right), namely d± \\ r\± || r2± (which we call "zzz") and 
d± _L ri± || r 2 ± ("zyy"). We have also measured the orientation-averaged quantity ("ave") defined 
as tfg ve (6;z±; \Ri±\, \R2±\) = / dR l± JdR 2 ±%:(e;z±;RiuR2i.), where / dR iX stands for integra- 
tion over the orientations of R[±. The lattice version of this equation is easily recovered for even 
(integer) values of the transverse sizes; in our particular case, |r,-j_| = 1, we have to use a sort of 
"smearing" procedure, averaging nearby loops as depicted in Fig. [T] (right). 

Since we are interested in the limit T — > °°, we have to perform it on the lattice by looking 
for a plateau of ^ plotted against the loop lengths L\^. On a 16 4 lattice it is difficult to have a 
sufficiently long loop while at the same time avoiding finite size effects, and at best we can push 
the calculation up to L = 8. Nevertheless, our data show that % is already quite stable against 
variations of the loop lengths at L\ 2 — 8 (at least for 6 not too close to 0° or 180°, where it is 
expected to diverge due to its relation with the static dd potential, see [0, ||) and so we can take the 
data for the largest loops available as a reasonable approximation of 

We have considered the values d = 0, 1,2 for the distance between the centers of the loops: as 
expected, the CFs vanish rapidly as d increases, thus making the calculation with our simple "brute 
force" approach very difficult at larger distances. 

4. Comparison with analytical results and the problem of total cross sections 

As already pointed out in the Introduction, numerical simulations of LGT can provide the 
Euclidean CF only for a finite set of 0-values, and so its analytic properties cannot be directly 
attained; nevertheless, they are first-principles calculations that give us (within the errors) the true 
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Figure 2: (Left) Comparison of the lattice data to the SVM expression (4. 1 ), with Kswm calculated according 
to [FTO] (solid line), or determined through a best-fit (for the "zzz" and "zxy" cases only, dotted line) at d = 1. 



(Right) Comparison of the lattice data to the ILM expression ( |4.2[ ) with ATilm calculated according to [|8|] 
(dotted, dashed, and solid line, corresponding to different values of the model parameters), or determined 
through a best-fit (sparse dotted line) for the "zzz" and "zyy" cases at d = 1. 

QCD expectation for this quantity. Approximate analytic calculations of this same CF have then 
to be compared with the lattice data, in order to test the goodness of the approximations involved. 
This can be done either by direct comparison, when a numerical prediction is available, or by 
fitting the lattice data with the functional form provided by a given model. The Euclidean CFs we 



are interested in have been evaluated in the Stochastic Vacuum Model (SVM) [jllj], in the Instanton 
Liquid Model (ILM) [|l2|, ^J, and using the AdS/CFT correspondence Jl3|]: the comparison of our 
data with these analytic calculations is not, generally speaking, fully satisfactory. 
In the SVM [O] the Wilson-loop CF is given by the expression 



(SVM) 



(0) 



-exp 



1 



-^SVMCOte 



1 (2 
+ 3 CXP V 3^ SVMCOt0 



1. 



(4.1) 



where A'svm is a function of z±, Ri± and 7?2_l only, given in []11|], that we have used to evaluate 
( fP| ) numerically in the relevant cases. The SVM prediction (44) agrees with our lattice data in 
a few cases, at least in the shape and in the order of magnitude, but, in general, it is far from 
being satisfactory, see Fig. ^ (left). The same conclusion is reached if one performs instead a one- 
parameter C^svm) best-fit with the given expression: the values of the chi-squared per degree of 
freedom (xlo f ) °^ tn * s an( ^ tne otrier fits triat we h ave performed are listed in Table [T[ 
The ILM predicts the following functional form of the CF Jl2[ g], 

KlLM 



(ILM) 



e) 



(4.2) 



sin0 ' 

in particular, a well-defined numerical prediction for .Ktlm has been obtained in [j8|]. The ILM 
prediction turns out to be more or less of the correct order of magnitude in the range of distances 
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Table 1: Chi-squared per degree of freedom for a best-fit with the indicated function. 



considered, at least around = § , but it does not match properly the lattice data; the same mismatch 
is seen also in a best-fit with Eq. (4.2), see Fig. || (right). Moreover, the ILM prediction seems to 



overestimate the correlation length which sets the scale for the rapid decrease of the CF with the 
distance between the loops: this is also supported by the comparison of the prediction for the 
instanton-induced dd potential Vdd with some preliminary numerical results on the lattice ||]. It 
is worth noting that largely improved best-fits (see Table |l|) are obtained by combining Eq. ([OJ) 
with the functional form < ^ pert ^(0) = A r pert (cot 0) 2 , corresponding to the leading-order result in 
perturbation theory 00, El, into the following expression, ^ (ILMp) (0) = §^ + £ I ' LMp (cot0) 2 . 

Finally, we have tried a best-fit with the expression obtained through AdS/CFT for the jV = 4 
SYM theory at large N c , large 't Hooft coupling and large distances between the loops []l3[|: 

^(AdS/CFT) ( Q \ = \ + Ki cot e + ^ 3 cos cot 1 - 1 . (4.3) 
[ sin J 

Taking into account that it is a three-parameter best-fit, even this one is not satisfactory: best-fits 
with QCD-inspired expressions with only two parameters, like, e.g., the ILMp expression [or some 
appropriate modification of the SVM expression (|Q|)] give smaller Xdof ( see Table |l|). 

As an important side remark, we note that our data show a clear signal of C-odd contributions 
in dd scattering, which are related through the crossing-symmetry relations [^] to the antisym- 
metric part of ^(0) with respect to = |. An asymmetry is present in the "zzz" and "zyy" 
tranverse configurations (^J ve is trivially symmetric), thus signalling the presence of odderon con- 
tributions to the dd scattering amplitudes. Although these C-odd contributions are averaged to zero 
in meson-meson scattering (at least in our model), they might play a non-trivial role in hadron- 
hadron processes in which baryons and antibaryons are also involved. 

As we have said in the Introduction, the main motivation in studying soft high-energy scatter- 
ing is that it can lead to a resolution of the total cross section puzzle. From this point of view, a 
satisfactory comparison of the lattice data with the SVM or the ILM would not have helped, since 
they yield constant or vanishing cross sections at high energy, as it can be seen by using Eqs. ( [Q| ) 
and the optical theorem. An ambitious question that one can ask at this point is if the lattice data are 
compatible with rising total cross sections. An answer can in principle be obtained by performing 
best-fits to the lattice data with more general functions, leading to a non-trivial dependence on en- 
ergy. This approach requires special care, because of the analytic continuation necessary to obtain 
the physical amplitude from the Euclidean CF: one has therefore to restrict the set of admissible 
fitting functions by imposing physical constraints (e.g., unitarity). 

In this framework, a possible strategy is suggested by the improvement of best-fits achieved 
with the ILMp expression: the idea is to combine known QCD results and variations thereof. 
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As an example, one could consider exponentiating two-gluon exchange and the one-instanton 
contribution (i.e., the ILMp expression), and supplementing it with a term which could yield a rising 
cross section, e.g., ^ nse ^ = expl^^^lexpjAnse^) (f — $) 4 (cot d) 2 } — 1. Such an expression 
yields an amplitude respecting unitarity if A r i se (&) > for large b = \z±\, and leads to the limit 
behaviour ~ (logs) 2 allowed by the Froissart bound if A^ib) ~ b~ 4 for large b. 



Another possible strategy is suggested by the AdS/CFT expression (4.3): one can try to adapt 



to the case of QCD the analytic expressions obtained in related models, such as jV = 4 SYM. 



As it has been shown in [15], by combining the knowledge of the various coefficient functions 
Ki in (4.3) at large b [13] with the unitarity constraint in the small-ft region, a non-trivial high- 
energy behaviour for the dd total cross section in JV = 4 SYM can emerge (including a pomeron- 
like behaviour a ~ s 1 / 3 ). Although of course Eq. (JO) is not expected to describe QCD, it is 



sensible to assume in this case a similar functional form (basically assuming the existence of the 
yet unknown gravity dual for QCD). Assuming moreover that the known power-law behaviour of 
the Kfs (expected for a conformal theory) goes over into an exponentially damped one (expected for 
a confining theory), in particular Kt, — > ce~ tlb , one obtains a rising total cross section proportional 
to the limit behaviour ~ (logs) 2 . 

It seems then worth investigating further the dependence of the CFs on the relative distance 
between the loops, as well as on the dipole sizes, as they could combine non-trivially with the 
dependence on the relative angle: these and other related issues (including the above-mentioned 
more general best-fits) are currently under study [|I6|], and will be addressed in future works. 
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